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^— ■>) ' We present straightforward proofs of estimates used in the adiabatic approximation. The gap 

, dependence is analyzed explicitly. We apply the result to interpolating Hamiltonians of interest in 

quantum computing. 
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I. INTRODUCTION 

m 
>. 

I The quantum adiabatic approximation has a long history, going back to Born and Fock [l| early 

in the development of quantum theory. Recently, the realization that the adiabatic approximation 
could be used as the basis for a method of quantum computation has generated a resurgence 
of interest in this topic. 
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Despite the existence of an extensive literature [3, 0, H @, on rigorous proofs of estimates needed 
to justify the adiabatic approximation, doubts have been raised about its validity 0] leading to 
' confusion about the precise conditions needed to use it [1, [13, |Tl|. In part, this is because some 

Oh. papers emphasize different aspects, such as the asymptotic expansion, the replacement of the 

requirement of a nondegenerate ground state by a spectral projection separated from the rest 
C I of the spectrum, dependence of first order estimates on the spectral gap, and even extensions 

to systems without a gap (T^ . Moreover, the qualitative gap conditions frequently presented in 
elementary texts [H, [lj| are known to be insufficient. 
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In this paper, we present a straightforward, yet rigorous, proof of the asymptotic estimates in a 
form which makes explicit the relevance of the gap. The proof given here is based on the one in 
Ref. S with the modifications introduced in Refs. |l3 and ilB- The main idea is to consider the 
physical time evolution as a perturbation of the adiabatic time evolution. This leads naturally to 
an integration by parts formula which is the main technical tool. 

Roughly speaking, the adiabatic approximation says that when a Hamiltonian changes slowly in 
time, the corresponding time evolution approximately preserves spectral subspaces. In particular 
an eigenstate tp{^) evolves with high probability to the eigenstate ilj{t) of the instantaneous Hamil- 
tonian H{t) if the energy curve E{t) does not come too close to any other energy level of H{t). To 
make this precise, one considers the family of Hamiltonians Hit/r), where t should be thought of 
as the "external" clock time and the parameter r defines a time scale which allows tuning of the 
speed of change of the Hamiltonian. It gives rise to a scaled time s = t/r. Using this notation 
time evolution is defined in terms of the family of time dependent Schrodinger equations 

idt4>r{t)=H{tlT)4>r{t) (1) 
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Figure 1: Energy band bordered by gaps. Shaded areas correspond to spectrum of H{s). The projection 
P{s) is associated with a{H{s)) n [fe_(s), 6+(s)], and Q{s) = 1 — P{s). 



or equivalently 

i9,Vr(s) =rff(s)^,(s). (2) 

In the most common scenario, the adiabatic theorem is an asymptotic expansion in 1/t for the 
error involved in estimating the time-evolved ground state of H(0) by the ground state of H{s). 



Remark: Allowed time dependence. The adiabatic theorems give results when the time dependence 
of the Hamiltonian in the unsealed time t is of the form H{t/T) for some r- independent family of 
operators H{-). 

It is important to understand that this excludes Hamiltonians with several independent time 
scales. In particular, the example by Marzlin and Sanders fS] which involves a Hamiltonian 
H{t/T,LL}()t) with fixed frequency uq and big time scale r does not belong to this class. To see 
this, notice that in terms of the dimensionless time s — t/r, the Hamiltonian is H{s,luots) which 
is T-dependent unless wot is taken to be constant. 



In adiabatic quantum computation [2, ll7| , the Hamiltonians of interest are usually interpolating 
Hamiltonians of the form 

H{s) = [l- f{s)]Ho + f{s)H,. (3) 

Generalizations of the form 

H{s) = [1 - f{s)]Ho + f{s)H, + k{s)H2 (4) 

with fc(0) — fc(l) = have also been considered In most applications of the adiabatic theorem 
to quantum computation, there is another parameter, n, which describes the size of the system, and 
one is interested in the behavior of a family of quantum systems when n is large. In particular, one 
is often interested how certain quantities, such as the running time of a computer program, grow 
(or decrease) with n. This brings a new element into the discussion of the adiabatic approximation 
and requires a careful analysis of the error bound in terms of the gap of a family of quantum 
systems. We present such estimates in Theorems [3] and |4] of this paper. 

The essential assumption in the adiabatic theorem as it is presented here is that the spectrum 
(7{H{s)) has a band associated with the spectral projection P{s), which is separated by a gap g{s) 
from the rest. 

This setting leads to an adiabatic approximation where the error terms are 0(l/r'^) for some q > 1. 
There is a weaker form of the adiabatic theorem, where one does not require a spectral gap due to 
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Avron and Elgart [ij] and Bornemann [l9|. In this case the estimate on the error term is o(l) as 
T goes to infinity. 

There are several different viewpoints for a discussion of the quantum adiabatic theorem. Each 
one offers interesting insight. We shall only mention them briefly, since the main objective of this 
note is to give a short and concise proof of the adiabatic theorem in a setting which seems natural. 

Berry [20|] pointed out that time evolution of a quantum system in the adiabatic limit becomes 
geometric. Simon plj noted that Berry's discovery can be interpreted as a parallel transport in a 
vector bundle with Berry's phase as the corresponding holonomy. Kato [F] used the adiabatic time 
evolution as a powerful tool in his perturbation theory of linear operators. Born and Oppenheimer 
[23 | used the adiabatic approximation in order to separate fast and slow motion in molecules 
thereby explaining the qualitative picture of their spectra. Niu and Thouless 23] and Avron and 
Seller J3| showed that conductance in quantum Hall systems is defined in the adiabatic limit and 
is a topological number. 



II. STATEMENT OF RESULTS 

In the following, we will use the following notational conventions: The letter P denotes an orthogo- 
nal projection and at the same time its range. [A, B] is the commutator AB — BA of the operators 
A and B. The adjoint of A is denoted as A'' . Furthermore, we shall use regularly the following 
well known facts. 

• Let r be a positively oriented loop in the complex plane circling the spectrum associated with 
an orthogonal eigenprojection P of a self-adjoint operator H. Then P ~ ~{2m)~^ §-^{H — 
z)~^dz. 

• Let P{s) be a smooth family of orthogonal projections and Q{s) = 1 — P{s) the projection 
on the orthogonal complement. Differentiating the relation P'^{s) = P{s) gives 

Pis) = P{s)P{s)Q{s) + Q(s)P(s)P(s). (5) 

An operator of this type is called "off-diagonal." 

A summary of notation is given at the end of the paper (Table [l|. 

Let / C M be an interval with G I and H{s), s G /, a family of Hamiltonians. We are interested 
in the time-dependent Schrodinger equations ([1]), where r > defines a time scale. A change of 
variables converts the family of equations ^ to the equivalent family Wc will denote the 
unitary time evolution associated with Eq. ([2]) by Ut-{s). 

Throughout the article, H is assumed to fulfill the following conditions. 

(i) H{s), s Cz I, are self-adjoint operators on a Hilbert space Ti, with an s-independent domain 
D. 

(ii) H{s) is a fc > 2 times continuously differentiable map from / to the space B{D, H) of bounded 
linear operators from D to Ti. equipped with the graph norm of H{0). 

(iii) H{s) has gaps in the spectrum, and P{s) is the spectral projection on a band bordered by 
gaps, i.e., there are two real- valued, continuous functions and and a number g > 
such that 

dist({6+(s), (s))) >g {s e I), 

and P{s) is associated with the nonempty band (t{H{s)) D [6„(s), 5+(s)], s £ I. 
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Conditions (i) and (ii) imply that the resolvent s i-^ {H{s) — z)~^ as a map from / to B{Ti) 
is k times continuously differentiable and ensure the existence and uniqueness of the propagator 
Ur{s) if A; > 1 (see, e.g., Ref. Ill): Ur{s) is unitary, strongly continuous in s, maps D onto D, 
Ur{0) = idw, and for each ip e D, ■0r(s) := Ur{s)ip S C^{I,H) and satisfies Eq. Condition 
(iii) is illustrated by Fig. 1. 

The first two adiabatic theorems give an estimate of the probability of finding a state "0 initially in 
P(0) outside P(s) at time s and compare Pt{s) := Ur{s)P{0)Ur{s)'' to P{s) for large time scales 
r. Pt(s) is the time-evolved projection, whereas P{s) is the instantaneous spectral projection 
associated with the band of the spectrum of interest. 

Let us remark that although bounds on the transition probability follow directly from bounds on 
— as can be seen from relations ([8]) and Q below, they are of sufficient importance 

that we state them separately in the theorems that follow. 

Theorem 1 (First adiabatic theorem). Assuming that the general conditions on H stated above 
are fulfilled, the following holds. For ij] € P(0) the transition probability is of order 0(l/r^), i.e., 

{Ur{s)i;, [1 - P(s)][/,(s)V)= 0{1/t^) 

and 

\\PAs)-Pis)\\=Oil/r). 
Both estimates hold uniformly on compact subsets of I . 

Theorem 2 (Second adiabatic theorem / switching theorem). If in addition to the general condi- 
tions we assume that H is compactly supported in an interval ]0, l[c /, the above estimates can be 
improved to 

(t/.(s)^, [1 - P{s)]UAs)^) - (^(l/r^e^-i)), 
\\PAs)-P{s)\\=0{1/t'^-'), 

uniformly in I\ ]0, 1[. 

Remarks concerning Theorem\^ If H is compactly supported in an interval ]a, a;[c /, the previous 
estimates change into the following: For S P{ci), 

{Ur{s; a)^, [1 - P{s)]Ur{s; a)^) = 0{1/t^( 
\\Uris;a)Pr{s)Ur{s;a)^ ~ P{s)\\ = 0{1/t^-^), 

uniformly in The miitary Ur{s]a) denotes the propagator of the Schrodinger equation 

for the inital value a. 

For even more regular Hamiltonians, different methods of proof can be used to show that the 
transition amplitude in Theorem [5] becomes exponentially small in r. For example, when the 
family of Hamiltonians H{s) is in a Gevrey class, one can apply methods of phase space tunneling 
[25j . or use Nenciu's technique [1^, exponentially small bounds are also given in Refs.iand0. 

Theorem 3. Suppose that the spectrum of H{s) restricted to P{s) consists of m{s) eigenvalues 
(each possibly degenerate, crossing permitted) separated by a gap g{s) from the rest of the spectrum 
of H{s), and H, H, and H are bounded operators. (This assumption is always fulfilled in finite- 
dimensional spaces.) Then for ijj G P(0) and any s G /, 
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and /lu.b. is a shorthand for /(O) + /(s). 



Theorem 4. Suppose that the assumptions of Theorem^ are satisfied and H is three times differ- 
entiable. Let h{s) := max(| |ij(s)| |, ||ff(s)||, \ \{d?' / ds^)H{s)\\) . Then there is a constant C (possibly 
depending on m) such that Eq. ^ holds with 
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A{s) = — 
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III. NOTATION AND IDENTITIES 



A. Adiabatic time evolution and wave operator 



The proof of the adiabatic theorems is most easily accompHshed with the aid of an ideahzed time 
evolution, the adiabatic time evolution introduced in Ref. mapping the spectral subspace PifS) 
onto P{s). This is an adaptation of an idea by Kato to the present situation [5]. 

From the gap condition (iii) and {H{-)-z)-^ G C'il, B{H)), it follows that P{s) is in C'=(/, B{n)) 
as well. Define the adiabatic Hamiltonian associated with P and the time scale r by 

H^{s) ■.= H{s) + -[P{s),P{s)], sel. 

T 

satisfies the conditions (i) and (ii) with differentiability of order fc — 1. Existence of the 
corresponding time evolution U^{-) is ensured for fc — 1 > 1. 

The adiabatic time evolution generated by tH^^s) is ideal in the following sense. 
Lemma 1. U^{-) satisfies the intertwining property 

U^{s)P{Q)=P{s)U^{s), sel. (7) 



Sketch of proof. The strategy of the proof is to show that both sides satisfy the same initial value 
problem X{0) = P{0), X{s) — —iTH^{s)X{s) using that P is off-diagonal, see identity ([5]) (see 
also Ref.i, p. 38). ■ 



In the following, we will omit the s-dependence and use the notation Q = 1 — P, Pq, Qo as a 
shorthand for Q{s) = 1 - P(s), P(0) and 0(0). 

We prove the theorems by comparing the real time evolution Ut{s) with the idealized time evolution 

U:^{s). The wave operator ^ris) := C/^(s) Ur{s) is a standard tool used in scattering theory to 
compare dynamics, here, Ur and U:^. We expect that adiabatic and real time evolutions are close, 
i.e., is close to the identity, for big r. The proof of the adiabatic theorems reduces to the 
estimation of the off-diagonal blocks of the wave operator since for tp G Pq, HV'H — 1; 



{Ur^,QUri^) = IIQC/^Po^lP = \\QU^ilrPo4'\? < WQo^rPoW^ 



(8) 
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and 

||P,-P|| = WUrPoUr^ -U^Po{U:^)^\\ = \\n,Po~Po^r\\ = HQo^rPo ~ Po^rQoW- (9) 

By a straightforward calculation fir satisfies the Volterra integral equation 

Qris) = 1 - / Kr{s')nr{s')ds' , (10) 

Jq 

where the kernel Kr is defined by 

Kris) ■.^U^is)^[Pis),Pis)]U^{s), 
and the integral here and in the following is the Riemann integral in the strong sense. 



B. Integration by parts lemma 



With a map X from / to B{T-C), we associate another map X[-] defined by 

X[s] U^is)^Xis)U^is). 

With this notation the integral kernel of the Volterra equation (fTUl) introduced in the last 
section is just [P, P][-] and the Volterra equation gives rise to integrals of the form J X[t]Y{t)dt. 
Such integrals can be integrated by parts. 



In the integral by parts formula (Lemma [3]), the twiddle operation plays a crucial role. It can be 
understood as a partial inverse of the add operation A i-^ [H, A] as follows. 

Lemma 2 (Twiddle operation). Let X be a map from I to B{Ti.) and r(s) a contour in C such 
that P(s) = -(27ri)-i/p^^^[i/(s) - z]-'^dz. Define 

X{s) := — / [H{s) - z]-^X{s)[H{s) - z]-^dz. 
27r» Jr{s) 

Then PXP = QXQ = and on D 

[H, X]=PX- XP 
P[H^, X]Q = PXQ, Q[H^, X]P = ~QXP. 

Proof. We drop the s-dependence. Let V be a contour in C lying outside the domain delimited by 
r and such that P and X can be written as integrals over V instead of F. Then 

= -tAi / - z)-^X{H - z)-i / [-^ - (if - z')-']dz'dz 

= -—^ i (H ^ z)-'X{H - z)-^ i —^{H - z)-^{H - z')-^dz'dz 
{2ttiY Jp 7p, z- z' 

"^—j j j^,{H - z)-^dzX{H ~ z')-^dz' 



{2m) 



— i P(H - z')-^X(H - zY^dz' = PX. 

2711 Jp, 
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This implies PXP = = QXQ. The commutator [H,X] can be computed on _D D Tl{X), 

[H, X]^ — (f[H ~ z,{H ~ z)-^X{H - z)-^]dz 

= — i{X(H - z)-^ -{H - z)-^X}dz = -XP + PX. 
2-Ki Jr 

The second pair of identities foUows with 

P], 1]Q - P]P)XQ - PX{Q[P, P]Q) = 0, 

and similarly Q[[P, P], X]P = 0. ■ 

Lemma 3 (Partial integration). Let X e C^{I,B(n)), Y e C"{I,B{H)). Then X defined as 
above is in C^{I,B{T-l)) and for s e /, 

r PoX[s']QoY{s')ds' = --(poX[s']Qoyis')^ - [ PoX[s']QoYis')ds' 
Jo TV Jo 

PoX[s']QoY{s')ds' 





A similar equality holds with Pq and Qq interchanged and an overall change of sign on the right- 
hand side. 

Remark. The notation X[s] is a shorthand for adjoining {U^)'' and [/^ after taking into account 
other symbols on X, e.g. X[s] = U^{s)'' X{s)U^{s). 

Proof. It suffices to notice that X[-] is strongly differentiable and 

^PoX[s]Qo = Po(ir[i/^,X][s] +i[s])Qo = iiTP[H^ , X]Q)[s] + PaX[s]Qa 
as 

= Po{iTX[s]+X[s])Qo, 
where the intertwining property has been used in the form 

iPZQ)[s] = P[s]Z[s]Q[s] = PoZ[s]Qo 
for Z : I ^ B{H). ■ 

IV. PROOF OF THEOREM [l] AND THEOREM [2] 

Using the intertwining property, we see that the integral kernel of the Volterra equation ITUl) can 
be written as 

Kr{s) = [PP - PP)[s\ = {QPP - PPQM - QoP[s\Pq - PoPWOo. 
As a consequence, the Volterra equation (jlOp written for Q^VlrPo becomes 



Qonr{s)Po = - / QoP[s']Ponr{s')Pods' . 

Jo 
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We can apply Lemma [3] and integrate by parts 

Qo^^r(s)Po = -^{QoPW]Po^lr{s')Po\o - ^ QoPW]Po^ris')Pods' 

QoP[s']PaPW]Qon,is')PQds'y (11) 

Therefore ||Qo^^t^o|| — 0{1/t). The same method gives the bound llPof^rQoH = 0(1/t), so that 
Theorem [jfoUows with Eqs. ^ and 

Now, suppose that H is compactly supported in ]0, 1[. The Volterra equation pUj) shows flris) = 1 
for s < and ilr{s) — ilri^) for s > 1; hence, it is enough to prove the statement for s = 1. In 
Eq. (|lip . for s > 1, the boundary terms vanish. The basic idea to prove Theorem [5] is to estimate 
the two remaining integrals by iterating the integration by parts. However, doing so naively 

for the second integral with X{s) = P{s) and Y{s) = P[s](5o^^r (s)^o leads to derivatives of 
P[s] = U^{s)'' P{s)U^{s) and thus to powers of r. This can be avoided by inserting the expression 
for QoVItPq recursively into the second integral. As a price to pay, we have to consider multiple 
integrals. 

Lemma 4. Let < m < k. Then for all j eN and Xi, ...,Xj e C^"(]0, 1[, B{l-L)) the integral 

Aj{Xi,...,Xj) ^Qa dsiXi[si]Qo....Qo dsjXj[sj]PQnrisj)Po 
Jo Jq 

is of order 0(l/r™). 

Remark. In this notation, (3of^r(l)-Po = ^^i{P)- 

Proof. By induction over m. The statement for to = is obvious since all Aj are bounded 
as functions of t. For the induction step TO<fc— 1— s-to + I assume that Xi, - .^Xj are in 
C™^^(]0, 1[,B{TC)). Suppose first that j > 2. We integrate by parts with respect to the variable 
Sj and use the vanishing of boundary terms, 

QoXj[sj]PQnr{sj)Podsj 



QoXj[Sj^l]PonriSj-i)Po - / QQXj[Sj]PoilriSj)PodSj 

Jo 

' QoX,[sj]PoP[s,]Qonr{sj)Podsj^ 

QoXj[Sj-i]PQ^r{S]-l)Po - / QoXj[Sj\PoVLr{Sj)Pods.j 

Jo 



+ r ' dsjQo{X,PP)[sj]Qo r dsj+iP[s,+i]Ponr{sj+i)PQ 

Jo Jo y 

Inserting this into the definition of Aj, we get 

Aj{Xi, ...,Xj) - ^Aj-i{Xi, ...,Xj-2,Xj-iQXj) - Aj{Xi, Xj-i, Xj) 

— Aj^i{Xi, Xj^i, XjPP, P) 



This remains valid for j = 1 if we set :— 0. The new variables Xj^iQXj, Xj, XjPP, and P 
are in C™(]0, 1[, B{Ti)) (here we use to < fc — 1). Hence, 

A,{X,,..,X,)='-o(^)=0 



rm+l 
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Proof of Theorem\^ Suppose that H is compactly supported in ]0, 1[. Then P belongs to 
Co"^(]0, Lemma H implies Qof^r(l)^o = ~M{P) = 0{1/t^-^). The lemma still holds 
with Q and P interchanged so Fo^^t(1)(3o = 0{1/t^~^). Together with Eqs. ^ and this 
concludes the proof. ■ 



V. DEPENDENCE ON THE GAP 



Up to now, we have examined the dependence of the error terms on the time scale r only. In order 
to give error bounds with explicit gap dependence, we need to know how P decorated with dots 
and tildes behaves. The aim of the following series of lemmas is to generalize the bounds given in 
Ref. [13, pp. 8-9, to the case where P projects on more than one eigenvalue. These allow to prove 
Theorems |3] and SI The reader interested in projections P{s) associated with a single (possibly 
degenerate) eigenvalue can find a simpler proof of the inequalities stated in Lemma |8l up to a 
constant factor, in Ref. l27l. 

We will assume throughout this section that all operators are bounded and that PHP has discrete 
spectrum. This is always fulfilled in finite-dimensional spaces. We write 

m(s) 

P{s) = ^ P,{s), H{s)P,{s) = \j{s)P,{s), 

where Pj{s) are the projections associated with the m(s) eigenvalues Ai(s), ... ,Am(s)(s). Note 
that each eigenvalue might be degenerate and the eigenvalues are allowed to cross. Furthermore, 
let 

g{s) := dist(^{Ai(s), A„(s)}, a{H{s))\{Xi{s), .., X,n[s)is)}^ 
denote the gap in s. In the following, we will drop the s dependence and use the notation 

= [His) - z]-\ = Q{[H{s) ~ z] \q}-'Q. 

Note that the reduced resolvent Rz is well defined even if z is one of the eigenvalues Xj and 

II^A.II < 1/5. 

Lemma 5. Let X^A^B he hounded operators. Then 

ni 

X^- Y^iPjXRx^ + Rx^XPj), (12) 

^ <j) RzARzBRzdz = (P - Q){AB + AB - AB). (13) 

Proof. By Lemma d [H,X] = [P,X]. This implies 

Pj[H, X]Q = X^PjXQ - PjXQHQ = PjXQ. 

Hence, 

PXQ = J2 P]XQ = - ^ PjXRx^ . 
j j 

Similarly, QXP = - R\,XPj. Since by Lemma [HX is off-diagonal, Eq. ^ follows. Now, let 

G{A,B) ■.= J-<{ RzARzBRzdz. 
2m Jr 
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Then 

[H, G(A, -S)] = ^ - ^' RzAR,BR,]dz ^ AB- AB. 
An argument similar to the reasoning above fixes the off-diagonal part of G{A^ B), 

PG{A,B)Q ^ P{AB ~ AB)Q, QG^A, B)P = -Q{AB - AB)P. (14) 

The diagonal blocks of G{A, B) can be computed, using that z i—f R^ 'is holomorphic in the domain 
enclosed by F and has i?^ its derivative with respect to z, 

gG(A, S)Q = V 7^ / -^—R,APjBR,dz + / R,AR,BR,dz 

= - H Rx,APjBRx^ + = -QiSQ. (15) 

PG(A, S)P = V / , ^ , ^ PjAR^BPkdz 

j.k ^ ■> 

+ 77- ^ T. 77T^ TT. ^dz RAPkBPi. 

^ 2m Jr (Xj - z){Xk - z){Xi - z) ' 



Similarly, 



Now, 



For j ^ fc, 



2TTi 7p {Xj — z) 



■^Rz — Rx 



Rz — —{Rxj — Rxk) — R\jR\k: 



2Tri Jy {Xj - z){Xk - z) Xj - Xk 

and /p[(Aj — z)(Xk — z){Xi — z)]^^dz — for any j, fc, I e {1, m(s)}, so that 

PGiA, B)P = J2 PjARx^Rx.BPk = PABP (16) 

Putting together Eqs. (ll4p - (fTB|) and making use of the off-diagonal character of any "twiddled" 
operator A, B, etc., gives Eq. ([Tc 



Lemma 6. We have P ~ H . Let X he a continuously differentiable map from I to BiT-L). Then 
X^X + {Q- P){PX + XP+ [H, X] - [P, X]). (17) 

In particular, 

QPP = Q{P + HP - PH)P. (18) 

Proof. First notice 
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Identity ([TT]) then follows with 

X = -^-^ / R,XR,Jz <p (RzXR^ - R.HR^XR, - R^XR,HR,)dz 



and Lemma m Since PX and XP are diagonal operators and [P, P] = 0, Eq. (fT8|) follows from 
Eq. nil) apphed to X ^ P. ■ 



Lemma 7. Let X : [0, 1] B{n). Then 



™||<y^M^. (19) 



Proof. Let e 

IIX0II = \\{PXQ + QXPM^ = WPXQc^W' + WQXPc^W' 

<max(||Plg||M|QlF||2 

Moreover, 

" ~ t 

Qxp = ~J2 R^^^Pj = (PX^Q) 

and for bounded operators \ \A\ \ ~ \ \A^\\, hence it is enough to prove 

\X\\ 



\\PXQ\\< 
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But this follows from \\PXQ\\ = \\{PXQ)'' {PXQ)\\^^^ and 

m 

WPXQW' - \\iPXQ)\PXQ)\\ = \\J2Rx,X^PjXRx,\ 

< mh\X\\\\X^\ ^ m^-^^. 
5 9 

Lemma 8. T/ie following series of bounds holds: 

\H\\ 



\P\\ < V^- 



9 9 

ii<3PPii<^m+M|nf 



5 5 5 5 

More generally, let TX :— X, DX :— X, then if H is (I + l)-times differentiable and h{s) :— 
maxi<j<;_|.i ||(^)-'^f(s)||, there exists a constant Ci, possibly depending on the numberm of eigen- 
values, such that 



y fc,Qi,...Q,>i: ]^<,<fc y 9 9 

Qi + ..+afc=i + l — 

i.e., every tilde pulls out a l/g and every dot an h/g. 
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The proof of the series of bounds is straightforward with Lemma IH] and Eq. ([7]). The estimate (PU)) 
is shown with an induction over I. 



Theorem[3]now is a simple consequence of Eq. pT|) and Lemma[71 If H is three times differentiable, 
iterating the integration by parts gives the expansion 



T T" 

1 r . r 1 ? ~~ . ~ . ~ (21) 

-— / Qo{PP)[s'W)P[0\P^~— Qo{PP + PP-PPP)[s']nris')Pods' 

''' JQ '''Jo 

+ \ f Qo{PP)[s] r {PP + PP)[s"]nr{s")Pods"ds'. 

Jo Jo 

Theorem |4] follows from this expansion and Lemma [H 

Comment on the traditional "adiabatic criterion. " How big must t be in order to ensure that the 
error HQo^^t^oH remains smaller than some constant e < 1? A frequently given answer is 

ll^(s)|| 

r > const sup — . 

o<s<i gysY 

This criterion comes from the first order term in the expansion (j2ip which we rewrite as 

||QoS2r(s)-Po|| < - ^(0) + —{s) + —. 

TVS a ) T 

However, C{H) in the second order term may depend on the gap. For a fixed Hamiltonian, one can 
always justify neglecting higher order times by taking r sufficiently large. However, in applications 
to quantum computation one is often interested in a family of Hamiltonians in which the gap may 
change. Indeed, in the situations considered in Theorems 3 and 4, the second order coefficient 
C{H) can be 0(1/5^). In such a case, this term is unbounded for vanishing g unless r = 0(1 /g^). 

There have been attempts toward a rigorous criterion based on an estimate of the form 



h 



1 h''+ 



\\Qonris)Po\\<C{— + sup ^^^) (22) 



u.b. 0<s'<s 9 



for (k + l)-timcs differentiable Hamiltonians (Ref. i27|, p. 2), where 

h{s) — max ||^^(s)||. This bound suggests that the transition amplitude can be kept small by 

choosing r of the order of g^+^Z'^. For large fc, this is very close to the traditional choice t ^ 
and much better than the scaling t ^ g^ obtained from Theorem 31 

The methods presented here do not allow to prove the bound ([22]) when k > 2 (compare to 
Theorem |4l). Trying to prove it by iterating the integration by parts fails; as we saw in Sec. IIV| a 
direct iteration of the integration by parts generates powers of r (see p. [5]). This was overlooked 
in Ref . |3 and later on corrected in Refs. [H and [H. However, it seems that a bound of type ((22)) 
can be proved using Nenciu's expansion, as sketched in Ref. [28|. 



VI. APPLICATION TO INTERPOLATING HAMILTONIANS 



In adiabatic quantum computation, the Hamiltonians of interest are interpolations between an 
initial Hamiltonian Hq with an easily computable ground state and a final Hamiltonian Hi whose 
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ground state encodes the solution to some problem. The running time of an "adiabatic algorithm" 
is closely related to the time r required to ensure a good agreement between real and adiabatic 
time evolutions and depends crucially on the minimal gap gmin- 

At first glance, Theorem [3] does not imply the general wisdom that r should be of order l/^^^jj^: 
if we estimate the integrals by the maximum of the integrand, we get a worst case dependence in 
l/ffmin- ^ closer look at the following example however shows that it is not a good idea to bound 
the integral in this way and that additional knowledge about the gap function q(s) does allow to 
extract the l/(?min behavior from Theorem[3l Moreover, as suggested in Refs. |23,l30, it is possible 
to improve the error dependence on the minimal gap by adapting the interpolation between Hq 
and Hi. 



Example: (See also Ref. ^29) Let |0), |1) denote an orthonormal basis of C^. For w G {0, 1}", let 
\w) = \wi)^---(^\wn),n:= (C2)®", 

Ho:=l-|6)(6|, Hi:=l-\u){ul 

where |6) — 2~"/^ ^^^^q j-j„ \w) and u S {0,1}". The element u is regarded as an unknown 
element to be searched for in a list of length 2" . We claim the following. 

(1) Let H{s) = (1 — s)Ho + sHi. Then g,nin = 2^"/^, and there is a constant C > such that 
||Qof^r(l)Po|| <C/(r5^J. 

(2) Let H{s) = (1 — /(s))ifo + f{s)Hi. It is possible to choose / : [0, 1] — > [0, 1] onto, monotone 
strictly increasing, such that ||(5of^r(l)^o|| < C/lT^min) for some constant C > 0. 

Proof. H{s) = (1 — s)Ho + sHi has a nondegenerate ground state separated from the first excited 
state by the gap 



g{s) = ^2-" + 4(l-2-")(s-i)2. (23) 



It is 1 at the boundaries s = 0, 1. We notice that for p > 1, g P{s)ds = 0{g]^^^) and ||ij|| = 
||i/i — HqW < 2. The first claim now follows from Theorem [3l 

Next, H{s) = [l-f{s)]Ho+f{s)Hi has gap g(/(s)) with g defined in Eq. Let / : [0, 1] [0, 1] 
be the solution of 

/(O) = 0, f{s) = kgP{fis)}, k - g-P{u)du. 



Then /(I) = 1. If p > 1, fc = 0(5^„^). Since 

H{s) = kgP{f{s))iHi ~ Ho), H{s) = k^g'^-' {f{s)) g(/(s)) {H^ - Ho), 
Theorem [3] gives 

\\Qonr{s)Po\\ < i(4fc + A:2^'[2.g2f-3(/(s)) \g{f{s))\ + 28g'P-%f{s))]ds 
k 

< 



-(4 + [2gP-^{u)\g{u)\+28gP'\u)]du 

(we used the change of variables u — /(s)). For p < 2, J g^^^ — 0{g^^^). The function g is 
decreasing on [0, 1/2] and increasing on [1/2, 1]; therefore for p < 2 

'9'-'\9\ = -^il-9'r^,!) = 0{9f-,^,). 
p z 







Hence, for 1 < p < 2, | |Qof^r(l)Po| I = 0{t-^ gl-.^g^^.^J = 0((r5„,i„)-i). This concludes the 
proof. ■ 
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More generally, the previous considerations concerning the gap dependence will apply provided 
that the gap function g{s) has the three following properties: 

- The gap at the boundaries {s — 0, s = 1) is of order 1. 

- The gap decreases strictly towards its minimal value gmin and then increases strictly. 

- If p > 1, then J^g-n-s)ds = 0(g^„P). 

The first two features are shared by most examples studied so far [13, HH . The third one needs to 
be checked more carefully, but it is reasonable to expect that J g~'P behaves better than (jj^ifn- 



Table I: Summary of notation. 


Symbol 


Meaning 


U-ris) 


real time evolution 




adiabatic time evolution 


X[s] 




Xis) 


(2«)-i §^^^^[H{s) - z]-'X[H{s) - z]-'dz 


nr(s) 


U:^{s)^ Ut{s) (wave operator) 




[H{s) — z]^'^ (resolvent) 




Q{[H{s) — z] \q}^^Q (reduced resolvent) 
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